A Nitsche-based cut finite element method for a fluid--structure interaction problem

We present a new composite mesh finite element method for fluid--structure interaction problems. The method is based on surrounding the structure by a boundary-fitted fluid mesh which is embedded into a fixed background fluid mesh. The embedding allows for an arbitrary overlap of the fluid meshes. The coupling between the embedded and background fluid meshes is enforced using a stabilized Nitsche formulation which allows us to establish stability and optimal order \emph{a priori} error estimates, see~\cite{MassingLarsonLoggEtAl2013}. We consider here a steady state fluid--structure interaction problem where a hyperelastic structure interacts with a viscous fluid modeled by the Stokes equations. We evaluate an iterative solution procedure based on splitting and present three-dimensional numerical examples.Comment: Revised version, 18 pages, 7 figures. Accepted for publication in CAMCo

Efficient white noise sampling and coupling for multilevel Monte Carlo with non-nested meshes

When solving stochastic partial differential equations (SPDEs) driven by additive spatial white noise, the efficient sampling of white noise realizations can be challenging. Here, we present a new sampling technique that can be used to efficiently compute white noise samples in a finite element method and multilevel Monte Carlo (MLMC) setting. The key idea is to exploit the finite element matrix assembly procedure and factorize each local mass matrix independently, hence avoiding the factorization of a large matrix. Moreover, in a MLMC framework, the white noise samples must be coupled between subsequent levels. We show how our technique can be used to enforce this coupling even in the case of non-nested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. We observe optimal convergence rates for the finite element solution of the elliptic SPDEs of interest in 2D and 3D and we show convergence of the sampled field covariances. In a MLMC setting, a good coupling is enforced and the telescoping sum is respected.Comment: 28 pages, 10 figure

A mixed finite element method for nearly incompressible multiple-network poroelasticity

In this paper, we present and analyze a new mixed finite element formulation of a general family of quasi-static multiple-network poroelasticity (MPET) equations. The MPET equations describe flow and deformation in an elastic porous medium that is permeated by multiple fluid networks of differing characteristics. As such, the MPET equations represent a generalization of Biot's equations, and numerical discretizations of the MPET equations face similar challenges. Here, we focus on the nearly incompressible case for which standard mixed finite element discretizations of the MPET equations perform poorly. Instead, we propose a new mixed finite element formulation based on introducing an additional total pressure variable. By presenting energy estimates for the continuous solutions and a priori error estimates for a family of compatible semi-discretizations, we show that this formulation is robust in the limits of incompressibility, vanishing storage coefficients, and vanishing transfer between networks. These theoretical results are corroborated by numerical experiments. Our primary interest in the MPET equations stems from the use of these equations in modelling interactions between biological fluids and tissues in physiological settings. So, we additionally present physiologically realistic numerical results for blood and tissue fluid flow interactions in the human brain